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Abstract 

Recent developments showed that light-cone parton distributions can be studied by investigating 
the large momentum limit of the hadronic matrix elements of spacelike correlators, which are known 
as quasi parton distributions. Like a light-cone parton distribution, a quasi parton distribution 
also contains ultraviolet divergences and therefore needs renormalization. The renormalization of 
non-local operators in general is not well understood. However, in the case of quasi quark distri¬ 
bution, the bilinear quark operator with a straight-line gauge link appears to be multiplicatively 
renormalizable by the quark wave function renormalization in the axial gauge. We hrst show that 
the renormalization of the self energy correction to the quasi quark distribution is equivalent to 
that of the heavy-light quark vector current in heavy quark effective theory at one-loop order. As¬ 
suming this equivalence at two-loop order, we then show that the multiplicative renormalizability 
of the quasi quark distribution is true at two-loop order. 
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I. INTRODUCTION 


Parton distribution functions (PDFs) characterize the structure of hadrons in terms of 
their fundamental constituents — quark and gluon partons. They play a crucial role in 
describing high energy scattering experiments involving hadrons. According to the QCD 
factorization theorem M, the hadronic cross section for simple processes such as inclu¬ 
sive deep inelastic scattering and Drell-Yan processes can be computed as a convolution of 
the partonic cross section and the PDFs. Whereas the former is perturbatively calculable, 
the latter are intrinsically non-perturbative and difficult to compute. In the held-theoretic 
language, the PDFs are dehned in terms of non-local light-cone correlators, which are intrin¬ 
sically Minkowskian and therefore forbid a direct lattice simulation. Currently the extraction 
of PDFs relies on a suitable parametrization and htting of the parameters to experimental 
lattice data liQ 


or 


Recent developments showed that the light-cone parton distribution can be di¬ 

rectly extracted from the large momentum limit of the hadronic matrix element of a spacelike 
correlator, which is known as the quasi parton distribution. The quasi parton distribution 
does not have a real time dependence, and thus can be simulated on the lattice. Moreover, it 
has the same infrared (IR) behavior as the light-cone parton distribution. Its connection to 
the light-cone parton distribution can be established as a perturbative matching condition, 
which can also be viewed as a factorization. In Ref. 13| , we derived the matching condition 


at one-loop level for the non-singlet quark distribution. We showed explicitly that the soft 
divergences are canceled both for the quasi and for the light-cone distribution. The collinear 
divergences are the same for both distributions as well. We also presented the one-loop 
matching factor that transforms the frame dependence of the quasi distribution into the 
renormalization scale dependence of the light-cone distribution. 

The factorization in Ref. [isj was given for the bare quasi quark distribution, where 
all helds and couplings entering the quasi distribution dehnition are bare ones. However, 
both the light-cone and the quasi distribution contain ultraviolet (UV) divergences and 
therefore need renormalization. In contrast to the renormalization of local operators, where 
all UV divergences can be removed by local counterterms, the renormalization of a non¬ 
local operator matrix element such as the parton distribution is rather distinct and less well 
understood. First, it contains UV divergences associated with coefficients that do not have 
a polynomial structure; second, there exist extra unphysical singularity structures, e.g. the 
light-cone singularities in the light-cone quark distribution or the axial singularities in the 
quasi quark distribution. These issues need to be properly addressed, in order to obtain a 
meaningful renormalized parton distribution. 

The renormalization properties of the light-cone quark distribution can be conveniently 
studied in the covariant Feynman gauge, as pointed out in Ref. j^. Nevertheless, the choice 
of light-cone gauge provides a simpler picture in that the distribution reduces to the biprod¬ 
uct of two quark helds, and the renormalization also becomes more straightforward. The 
disadvantage of this gauge choice is that the treatment of light-cone singularities is not clear 
a priori. There have been a number of papers discussing the regularization of light-cone 
singularities, or of similar singularities in more general axial gauges . Given the similar¬ 
ity between the quasi and the light-cone distributions, it is natural that the renormalization 
of quasi quark distributions also becomes more straightforward in the axial gauge. There 
is, however, an important difference between the light-cone and the quasi distribution with 


respect to their UV behavior. In Ref. [13|, we have pointed out this difference: For the 
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light-cone distribution, the momentum fraction is dehned as x = k~^/p'^, where k~^ and 
are the plus-momenta for the quark in the loop and the initial quark, respectively. The 
light-cone momentum fraction is restricted to [0,1], and the momentum {k~^,k~,k±) can 
roughly be power-counted as ~ (p"*", A^/p+, A) with A being a UV cutoff. The UV diver¬ 
gences in the light-cone parton distribution are then genuine UV divergences to which the 
usual renormalization procedure can be applied. In contrast, the momentum fraction for the 
quasi distribution is dehned a.s x = n ■ k/n ■ p, where n is a space-like vector, and is chosen 
to be along the 2 ;-direction in Ref. [l3[. The self energy correction to the quasi distribution 
has the same UV behavior as in the light-cone distribution, since all components of the loop 
momentum are integrated over. However, the vertex correction behaves differently, since the 
quasi momentum fraction x, or equivalently, the 2 ;-component of the loop momentum k^, is 
left unintegrated in the vertex correction. As discussed in Ref. [l3[, the quasi momentum 
fraction is no longer restricted to [0,1], it can extend from [—cx),-1-cxo]. Therefore leaving 
X unintegrated reduces the power of UV divergences, and leads to a UV convergent vertex 
correction at one-loop. Of course, the UV divergences appear when going beyond one-loop. 
However, as we will see, they can be sub divergences coming from subdiagrams only, not 
overall divergences. Consequently, all UV divergences in the (axial gauge) vertex correction 
can be removed by counterterms for subdiagrams from the interaction and therefore do not 
affect the renormalization of the quasi quark distribution; the renormalization of the quasi 
quark distribution then reduces to the renormalization of two quark helds in the axial gauge, 
provided that the non-trivial complication due to the choice of axial gauge does not matter. 
We explicitly show that this is indeed the case at two-loop order. 

We will focus in particular on the unpolarized non-singlet quasi quark distribution. By 
investigating the one-loop corrections to the quasi quark distribution, we hnd an equivalence 
between the self energy correction of the quasi quark distribution and the correction to 
the heavy-light quark vector current in heavy quark effective theory (HQET), so that the 
UV divergences in the former can be renormalized as the renormalization of the latter. 
Based on this observation and the one-to-one correspondence between the two-loop diagrams, 
we then assume an equivalence between the UV divergences of the two-loop self energy 
in the quasi quark distribution and of the two-loop corrections of the heavy-light quark 
current, and discuss the two-loop UV divergence structure of the vertex correction and the 
renormalization of the quasi quark distribution. To simplify our calculation, we choose 
dimensional regularization for UV divergences throughout this paper. In this way, the linear 
divergence present in a cutoff regularization is ignored. The result presented in this paper 
therefore shall be viewed as a hrst step towards a full understanding of the renormalization 
property of the quasi distribution on the lattice. 

The rest of this paper is organized as follows. In Sec. 2, we briefly review the one- 


loop calculation performed in the axial gauge in Ref. [13|, and discuss its connection to 


the Feynman gauge computation. We also present a renormalization formula for the quasi 
quark distribution, and give the one-loop renormalization factor that renders the quasi quark 
distribution UV hnite. In Sec. 3, we consider the two-loop correction and renormalization of 
the quasi quark distribution. Also some discussions on its renormalization beyond two-loop 
are given. Sec. 4 is our conclusion. 
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II. ONE-LOOP CORRECTION AND RENORMALIZATION FOR THE QUASI 
QUARK DISTRIBUTION 


Let us start by recalling the definition of the quasi quark distribution. For the unpolarized 
quark density, it is given as 


-ig dz'A^{0,0^,z')]ij{0)\P) , (1) 


r°° dz — 

q{x,fi,P^) = j —e*^^^(P|V^(0,0x,^)7^exp 

where we have chosen dimensional regularization for potential UV divergences, and fi denotes 
the renormalization scale dependence. 

In Ref. 13|], we computed the one-loop correction to the unpolarized quasi quark distribu¬ 
tion in the axial gauge = 0. The advantage of this gauge is that the gauge link in Eq. ([T]) 
reduces to unity, and the contributing Feynman diagrams become simple. The disadvantage 
is that it is non-covariant, and the prescription for the axial gauge singularity is not clear 
a priori. Since the dehnition Eq. ([1]) is gauge invariant, we can carry out the computation 
in any gauge. A convenient choice for the investigation of its renormalization properties is 
the covariant Feynman gauge, as mentioned in the Introduction. In the Feynman gauge, 
the contributing diagrams are given by those in Fig. [T] Actually there is a straightforward 
correspondence between the diagrams in Fig. [T]and the individual terms in the axial gauge 
result. For example, the diagrams in the second row of Fig. [1] correspond to the contribu¬ 
tion to the vertex correction from the term, the —{n^{p — ky + rp{p — ky)/n ■ {p — k) 
term, and the n‘^{p — k)^{p — ky/{{n ■ {p — k)y term of the axial gauge gluon numerator, 
respectively. 

Let us write down the contribution of each Feynman gauge diagram. We start from the 
vertex corrections, i.e. the diagrams in the second row of Fig. [T] The hrst diagram gives 

r n j 


Til = 


k^ 


. (2vr)^ 

= -ig‘^C, 


m 


m 


d^k 


u{p) 


2y 


+ 


{p- 

8mk^ 


ky 


u{p)6{x - —) 


(27r)^ i{}A — m?‘){p — ky {IP — rri^yij) — ky 


k^ 


u 


(2) 


where the quark mass m is introduced as in Ref. [13[ to regularize the collinear divergences. 

The above integrals can be computed in the same way as in Ref. [lH. After a Feynman 
parametrization and integration over k^ and /c^, we have the following contribution to the 
one-loop quasi distribution 
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FIG. 1: One-loop diagrams for quasi quark distribution in Feynman gauge, the conjugate diagrams 
are not shown. 


and leads to the one-loop correction (together with the conjugate diagram) 


qi2 = 


asCp 

271 


lx 1 X—l _ 1 
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1—X _ - 

—^ In ^ ^ 

1 — X X—l 1—X ’ 
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X > I , 

0 < X < 1 , 
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(5) 


Note that there is no UV divergence in the above results. This is because the momentum 
fraction x is left unintegrated. As x can extend between [—oo, +oo], leaving it unintegrated 
reduces the power of UV divergence, and therefore leads to a UV convergent result in the 
one-loop vertex correction. If one integrates over the momentum fraction x, as is done for 
the self energy diagrams in the hrst row of Fig. [H one will hnd a logarithmic UV divergence, 
in accordance with the usual UV power counting. 

The third diagram gives (note ri^ = —1) 
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d^k 

Wi' 
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which leads to the following result 


u 
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1—X 


, X < 0 


in dimensional regularization. 

Now let us look at the self energy diagrams in the hrst row of Fig. [H where we integrate 
over all components of the loop momentum. Since we are interested in the renormalization 
properties of the quasi distribution, we need the UV divergent part of the diagrams only. The 
hrst diagram is the usual Feynman gauge quark self energy diagram, it leads to the following 
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contribution to the wave function renormalization factor Zp'' in dimensional regularizat 
{D = 4 — 2e) and MS scheme 


ion 
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( 8 ) 


where S'^ = (47r)'^/r(l — e) accounts for the constant associated with 1/e in the MS scheme. 
The second diagram yields 
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d^k 
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In dimensional regularization, this leads to the following contribution to 

usCpS^ 1 
271 e 


Z^•2 — 


(9) 


( 10 ) 


To understand how the Wilson line self energy diagram contributes, let us expand the quasi 
distribution dehnition in coordinate space Eq. ([1]) to 0{g^). For convenience, let us separate 
the gauge link as 


So dz'n-A{0fi^,z') _ I'-pg-igT'” d2'n'A(0,Ox,^')jtpg-*s/o°° dz'n-A{0fi^,z') 


( 11 ) 


and associate the hrst factor on the r.h.s. to 'ijj{z) and the second to '0(0) order to form 
gauge invariant quark helds. It is easy to see that the Wilson line self energy diagram arises 
when the 0{g‘^) term comes from one of the path-ordered exponentials on the r.h.s. of 
Eq. (HTI) only. If each of the exponentials contributes an 0{g) term, this leads to the last 
diagram in Fig. [H Suppose we expand the second exponential to O^g"^), we have 


/ ^e‘”'-{p\4,(0,0^,z)Y-SM 


dz'dz"V{n ■ A(0, 0j_, z')n ■ ^4(0, 0j_, z"))0(O) \p) 
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ig‘^CF f d^k 1 

2 J (27ry {p — ky[n ■ {p — k)y 


( 12 ) 


The expansion of the hrst exponential yields the same contribution. This derivation also 
gives a natural prescription for the double pole at n- (p — /c) = 0. From the z', z" integration 
above, one can see that to achieve well-dehned Wilson line propagators in momentum space, 
one must have a +ie prescription for one 1/n- {p — k) propagator, and a —ie for the other. 
This leads to the prescription 1/[{n- {p — k)y + for the double pole l/(n • (p — /c))^. The 
above derivation is also formally consistent with the intuitive understanding of the third 
diagram in the hrst row of Fig. [T], where the Wilson line propagator 1/n ■ {p — k) and the 
corresponding self energy correction are similar to those for a quark matter held in the on- 
shell and zero momentum limit, its contribution can therefore be extracted analogously by 
taking the residue of the self energy correction in the zero momentum limit, which is given 
by 
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From Eq. ffT^ or Eq. fflSj) and the prescription for the donble pole, we obtain the following 
contribution to Zp'^ from the Wilson line self energy diagram 


Zl3 — 


o-sC 1 
271 e 


(14) 


Now we can write down the renormalization of quasi quark distribution. The quasi 
distribution can be renormalized as follows 

<iR(x,p‘) = [ ^Z(-)q(y,p‘), (15) 

J \y\ y 


where the subscript R denotes the renormalized quasi distribution, Z{x/y) is the renormal¬ 
ization factor required to cancel the UV divergence. The above relation can be rewritten 
as 

Qr{x,p") = j ^^[ZFZ{^)][Zp^q{y,p^)], (16) 

where we separate the held wave function renormalization constant Zp such that the factor 
in the second bracket is the quasi distribution with renormalized helds and can be computed 
using the standard Feynman rules, with the counterterms from the Lagrangian also being 
taken into account. The factor in the hrst bracket then cancels the remained UV divergence 
in the parton distribution. 

From the one-loop results above, we can see that there is no UV divergence in the vertex 
correction (the second row of Fig. [1]). The UV divergence arises only from the diagrams in 
the hrst row of Fig. [H the sum of which corresponds to the quark self energy in the axial 
gauge. Therefore, at one-loop level one needs to renormalize the self energy only, and the 
renormalization constant is given by 

Z(p) = S(p - 1)(1 + Z'D) = S(p - 1) - 1 

dvr e 

An interesting observation is that the diagrams in the hrst row of Fig. [T]are in one-to-one 
correspondence with the one-loop correction to the heavy-light quark vector current q'j^Q 
in the HQET. Also the Wilson line propagator 1/n ■ k is the same as 1/v ■ k in the heavy 
quark propagator, except that the vector n in the former is space-like, whereas the velocity 
V in the latter is time-like. However, as we see from Eqs. (^, (E]), flTOj) . ([12]), flT^ . ffTT|) . 
the UV divergences are actually identical for the one-loop diagrams in both cases (up to a 
factor of 2 since the contributing diagrams for quasi quark distribution are twice of those for 
heavy-light quark vector current). Therefore, the one-loop renormalization for quasi quark 
distribution requires the renormalization of self energy diagrams only, which is equivalent 
to the renormalization of the heavy-light quark vector current in HQET. 


III. RENORMALIZATION AT TWO-LOOP 

In this section we consider the two-loop correction to the quasi quark distribution. Fol¬ 
lowing our discussion of one-loop diagrams, we can divide the two-loop Feynman gauge 
diagrams into two classes: one corresponds, as the hrst row of Fig. [H to the two-loop dia¬ 
grams for the heavy-light quark current in HQET , or to the two-loop quark self energy 
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FIG. 2: Two-loop self energy diagrams for quasi quark distribution in Feynman gauge. The disjoint 
and conjugate diagrams are not shown. The two-loop Wilson line self energy diagrams are the same 
as the first four diagrams for quark self energy, but with gluons attached to the Wilson line. 


in the axial gauge; the other corresponds, as the second row of Fig. [H to the two-loop vertex 
correction. The hrst class diagrams are depicted in Fig. |2l where we do not show the disjoint 
diagrams with two one-loop diagrams separated from each other. Also conjugate diagrams 
are not shown. The Wilson line self energy diagrams are the same as the hrst four quark 
self energy diagrams, but with gluons attached to the Wilson line. 

As in the one-loop case, the two-loop diagrams in Fig. [2] are in one-to-one correspondence 
with the two-loop diagrams for the heavy-light quark current in HQET, we can therefore 
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FIG. 3: Feynman gauge two-loop vertex diagrams that contain sub divergences, the conjugate 
diagrams are not shown (to be continued). 


renormalize the UV divergences in Fig. [2] in the same way as is done for the two-loop 
renormalization of heavy-light quark current 22j. Based on the equivalence of their one- 


loop corrections, in the following we will assume they have the same UV divergences also 
at two-loop level, and focus on the two-loop renormalization property of the quasi quark 
distribution. The UV divergent two-loop contribution to the heavy-light quark current has 
been analyzed in Refs. 22, 2^. 

We now consider the two-loop vertex diagrams. At two-loop, the overall UV divergence 
of a vertex diagram comes from the region where both loop momenta, denoted as ki and 
^ 2 , are large. When either ki or ^2 is large while the other is small, a subdivergence arises. 
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FIG. 4: Feynman gauge two-loop vertex diagrams that contain sub divergences (continued) 


According to our definition of quasi distribution, the momentum fraction x = k^/p^, or 
equivalently fcf, is left unintegrated in the two-loop vertex correction. Since kf can go to 
±oo, leaving it unintegrated reduces the power of UV divergence, a simple power counting 
tells us that no overall UV divergence appears in the two-loop vertex diagrams. However, 
there exist sub divergences. In Figs. [3] and |4] we show all two-loop vertex diagrams that 
contain subdivergences. We are interested in the renormalization property of the diagrams, 
and therefore will focus on their UV behavior. Given the one-loop results in the previous 
section, the sub divergences in these diagrams can be computed straightforwardly. In Table U 
we list the UV divergent parts of all the two-loop vertex diagrams. 

Summing over all the contributions in Table HI we have 

p( 2 ) ^ + gcp - 2nf)q^^\ (18) 

izvre 

where we have used Tf = 1/2 for the SU{N) gauge group, Cp = (iV^ — l)/2iV, Ca = N, Uf 
is the number of active quark flavors, and = qu -|- qi 2 + ^ 13 . Although each contribution 
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UV divergence 

Fig. in 

UV divergence 

1, 2, 8, 9 (-Fconj.) 
3, 4, 7, 10 (-|-conj.) 
5, 6 (-l-conj.) 

11, 12 (-|-conj.) 

13 (-l-conj.) 

14 (-l-conj.) 

15 (-l-conj.) 

osCf Al) 

47re y 

^^{CF + CAmii + m) 
-^(2gii+gi2) 
^(2911+912) 

0 

^sCf^ 

2-Ke 

27re *<13 

1 (-kconj.) 

2 (-kconj.) 

3, 4, 5 (-Fconj.) 

6 (-kconj.) 

7 (-kconj.) 

8 (-kconj.) 

9 (-l-conj.) 

10 (-kconj.), 11, 12 

- hCA)m 

=^i^(9n+9l2) 

^sCf^ 

27re yiS 

0 

S(|C4 - 


TABLE I: Results for the UV divergent parts of two-loop vertex diagrams of quasi quark distribu¬ 
tion. (-|-conj.) means including the contribution of conjugate diagrams. 


in Table |T] has a distinct divergence structure, the sum of them is proportional to the one- 
loop vertex correction This indicates that all UV divergences in the axial gauge vertex 
correction can be removed by counterterms from the interaction (for the UV counterterms in 
the axial gauge see e.g. Ref. jl^) and therefore do not affect the renormalization of the quasi 
quark distribution. The renormalization of the quasi quark distribution is then equivalent 
to the renormalization of the two separate quark fields in the axial gauge. Expanding the 
renormalization formula Eq. fllbp to two-loop level, we can write the two-loop renormalization 
factor for the quasi quark distribution as 

[(S)Vd(^ + j) + (Z“')T(»-l). (19) 

where a, b come from the renormalization of two-loop self energy diagrams in Fig. [21 As 
is well-known, the double pole term a/e^ can be derived from the single pole at one-loop. 
Suppose we write a renormalization factor Z as a series expansion 


Z 


i=l j=l 




( 20 ) 


For the one-loop self ener gy diagrams in Fig. [H the renormalization factor in a general 
covariant gauge is given as |23| is a gauge parameter, and the Feynman gauge is obtained 
for ^ = 1) 

Zi,i = = {-^Cf + (3 - OCf) + 2^Cf = 3Cf, (21) 


which is twice the renormalization factor of heavy-light quark current in HQET, since the 
contributing diagrams to Zi^i are twice those to the heavy-light quark current. 

The renormalization group equation then leads to the following relation j2ll^23 


1 22 4 9 11 

^2,2 = —^(—2Zi^i ~ ~^CaCf — TIjCf, 


( 22 ) 


where we have used the fact that at one-loop Zi i is independent of the gauge parameter 
The coefficient a in Eq. flT^ is therefore given by 


a — -{-Cp + —CaCf — u^Cf')- 


(23) 
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The single pole coefficient b can be read off from the two-loop result for the heavy-light 


quark current in HQET in Refs. [22, \2M, 125 


b = 


127 28 


10 
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-TT H- Uf. 

27 9 ^ 


(24) 


Beyond two-loop, it is still true that the vertex contribution contains no overall divergence, 
but sub divergences only, as can be seen from the power counting above. In the axial gauge 
where the gauge link becomes unity, all these sub divergences can be removed by counterterms 
from the interaction, provided that the non-trivial complication due to the choice of axial 
gauge does not matter. The renormalization of the quasi quark distribution then reduces to 
the renormalization of two separate quark helds in the axial gauge. This remains, however, 
to be checked by explicit computations beyond two-loop order. 

In the above discussion, we choose dimensional regularization for UV divergences. It 
is interesting to investigate the renormalization property of the quasi distribution in the 
presence of a cutoff regulator used in Ref. 1^ to mimic lattice setting. At one-loop, the 
cutoff regulator, denoted as A, introduces a linear divergence in gis and Z 13 [l3[. As discussed 
in Sec. 2 , the prescription for the axial gauge double pole at x = 1 is given by 1 /((I—x)^+£^). 
This leads to the following result for Z 13 


~ asCp,, .2 

Zi3 = —— InA- 

271 e 


(25) 


Since the vertex correction also contains a linear divergence, the one-loop renormalization 
factor for the quasi distribution now becomes 


Z(,) = 4(„ - 1)(1 - ^(31n - ~)) 

Att e 


asCp 


A 


27r (1 — rjy -f ’ 


(26) 


where the e-dependence drops out when integrating over a smooth function of p, as discussed 


in Ref. [13|. We do not discuss the renormalization with a cutoff regulator beyond one-loop. 


since the use of a cutoff regulator is ambiguous at higher-loop order. Moreover, there is no 
unambiguous way to implement a cutoff gauge invariantly in the continuum, although this 
can be achieved on a discretized lattice. 


IV. CONCLUSION 

We have discussed the renormalization property of the quasi parton distribution, focusing 
in particular on the unpolarized quasi quark distribution in the non-singlet case. In this case, 
the bilinear quark operator with a straight-line gauge link appears to be multiplicatively 
renormalizable by the quark wave function renormalization constant in the axial gauge. We 
explicitly showed that this is true at two-loop order. In the covariant Feynman gauge the self 
energy diagrams are in one-to-one correspondence with the loop diagrams for a heavy-light 
quark vector current in HQET at the same order, thus the renormalization of self energy 
correction can be carried out as renormalization of the heavy-light quark vector current 
in HQET. The vertex correction of quasi quark distribution contains sub divergences only, 
which can be removed by counterterm for subdiagrams. Such features are expected to hold 
also beyond two-loop order, but remain to be checked by explicit computations. 
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